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. Abstract 



Tests of the Standard Model and its hypothetical extensions require precise the- 
oretical predictions for processes involving massive, unstable particles. It is well- 
known that ordinary weak-coupling perturbation theory breaks down due to inter- 
mediate singular propagators. Various pragmatic approaches have been developed 
to deal with this difficulty. In this paper we construct an effective field theory 
for resonant processes utilizing the hierarchy of scales between the mass of the 
unstable particle, M, and its width, P. The effective theory allows calculations 
to be systematically arranged into a series in and T/M, and preserves gauge 
invariance in every step. We demonstrate the applicability of this method by 
calculating explicitly the inclusive line shape of a scalar resonance in an abelian 
gauge- Yukawa model at next-to- leading order in P/M and the weak couplings. 
We also discuss the extension to next-to-next-to-leading order and compute an 
interesting subset of these corrections. 



1 Introduction 



An important part of recent and future high-energy physics experiments belongs to the 
detailed investigation of the production of heavy unstable particles, such as W and Z 
bosons, top quarks, Higgs bosons or perhaps other new particles. Prom single resonant 
or pair production the masses and couplings of the particles can be determined with high 
precision, provided theoretical calculations are equally precise. Since the decay widths of 
the particles are often non- negligible, this includes a consistent treatment of finite-width 
effects beyond the narrow-width approximation. 

Several higher-order calculations involving unstable particles have been performed in 
recent years, in particular for the line shape of the Z boson [1], W^-pair production [2], and 
tt-production [3]. In these calculations the finite width of the particles has been treated in 
a variety of often pragmatic approaches. While this may be adequate for the present, it 
is certainly desirable to formulate a theoretical framework that would allow for systematic 
improvements of the accuracy of such calculations. Moreover, future precision experiments 
require that such a framework be developed. 

The difficulty with unstable particles is that they cause singularities in propagators, 
if the scattering amplitudes are constructed according to the rules of weak-coupling per- 
turbation theory. A well-known remedy of the singularity is resummation of self-energy 
corrections to the propagator, which results in the substitution 



The self-energy has an imaginary part of order M'^g^ ~ MF, where F is the on-shell decay 
width of the resonance, rendering the propagator large but finite. "Dyson resummation" 
sums a subset of terms of order (gr^M^/[p^ — M^-|-iMF])" ~ 1 (near resonance) to all orders 
in the expansion in the coupling g^. This procedure raises the question of how to identify 
all terms (and only these terms) required to achieve a specified accuracy in and F/M. 
The failure to address this question may lead to a lack of gauge invariance and unitarity 
of the resummed amplitude, since these properties are guaranteed only order-by-order in 
perturbation theory, and for the exact amplitude. 

Many of the current approaches to unstable particles take the restoration of gauge 
invariance as their starting point. One example is the fermion-loop scheme [4], which is 
based on the observation that the dominant contribution to the width of the W and Z 
gauge bosons comes from fermion loops. The prescription reads to include the fermion-loop 
corrections in propagators and vertices, so that gauge invariance is maintained, since all 
fixed-order g^N ^ terms (where N j is the number of fermion flavours) are included. Besides 
the restriction to gauge bosons the disadvantage of this scheme is that one-loop vertices 
must be computed even for a leading-order approximation. Another scheme [5] constructs 
a gauge-invariant non-local effective action, which can be matched onto the two-point 
functions of the underlying theory. The gauge Ward identities are satisfied by construction 
in this scheme. Both schemes have been implemented in four-fermion production mediated 
by H^-pair production solving not only the gauge-invariance problem but also capturing 
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some sophisticated features such as running of the couphngs. However, it is not evident 
how to extend them to a systematic approximation of the scattering amphtude in powers 
of and T/M. 

Other approaches exploit the presence of two different momentum scales in the pro- 
duction and decay of weakly interacting unstable particles: M, the mass of the unstable 
particle, and its width, F <S M, both related to the location of the resonance pole. In the 
pole scheme [6] the scattering amplitude is approximated by an expansion in F/M around 
the poles in the complex plane. The coefficients of this expansion must be gauge-invariant 
in every order. In pair production of unstable particles this approximation is referred to as 
"double-pole approximation" [2]. It has been used in the complete next-to-leading order 
(NLO) calculation for VF-pair production [2], where in one- loop radiative corrections the 
leading term in the F/M expansion is sufficient. 

The pole approximation can be considered as the first step towards a systematic ap- 
proximation scheme to the scattering amplitude based on the separation of scales. Within 
the (double) pole approximation the process naturally consists of production and decay 
subprocesses, connected by the intermediate resonance(s). One-loop corrections separate 
into two subsets: factorizable corrections to the hard production/decay subprocesses, and 
non-factorizable corrections, accounting for interference between different hard subpro- 
cesses. This separation does not occur at the level of individual Feynman diagrams. For 
example, the self-energy correction to the unstable particle contains both, factorizable and 
non-factorizable pieces, and the two have to be separated carefully to avoid double count- 
ing. This is expected to be complicated for multi-loop calculations, or for a calculation 
of F/M suppressed contributions, which lie beyond the pole approximation. One- loop 
non-factorizable corrections have been extensively studied in the double-pole approxima- 
tion. A few important theorems have been proved [7] , and explicit calculations have been 
performed [8] . We also mention an approach that constructs the expansion in of the res- 
onant propagator squared in the distribution sense [9] , also addressing the issues of gauge 
invariance and systematic expansions. 

In [10] it was suggested to use the scale hierarchy F <C M for constructing an effective 
field theory, from which hard modes with momenta of order M or larger are removed. 
The effect of the hard modes is included into the coefficients of the effective Lagrangian 
and corresponds to the factorizable corrections. Non-factorizable effects are reproduced 
by the dynamical modes of the effective theory. This provides a more precise definition 
of "factorizable" and "non-factorizable" that generalizes to higher loop orders. In [10] it 
was shown that this approach is equivalent to the double-pole approximation within the 
accuracy of the latter (i.e. one-loop, leading order in F/M). However, the idea should work 
beyond these approximations. 

In this paper we pursue this idea and develop the effective field theory approach to 
unstable particle production from a systematic point of view. We take the attitude that 
the scale hierarchy F <^ M is the characteristic feature of the process, and that all other 
issues such as resummation and gauge invariance will be an automatic consequence of any 
theory that formulates the rules for an expansion in g"^ and F/M correctly. We identify 
the factorization properties of the process to any order in F/M, and define the effective 
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production and decay vertices and the effective Lagrangian in dimensional regularization. 
Tfie matcliing corrections are computed by expansion of Feynman integrals in momentum 
regions [11]. In this way we identify all contributions to the scattering amplitude at a given 
order in and F/M (and only these). We believe that this approach solves the conceptual 
difficulties that have so far been associated with the perturbative treatment of unstable 
particles in principle; however, for complex scattering processes, the implementation in 
higher orders than NLO still requires difficult calculations. 

Here we shall describe our approach for the case of inclusive production of a single 
charged scalar resonance in fermion-fermion scattering in an abelian gauge model. The 
specific set-up of the toy model, and the restriction to inclusive scattering, imply that 
a minimal set of effective fields needs to be introduced, and only the forward scattering 
amplitude needs to be considered. This allows us to concentrate on the essential features of 
the approach, such as the identification of momentum modes, factorization of the scattering 
amplitude, construction of the effective Lagrangian and matching, and the calculation of 
the scattering amplitude in the effective theory. We compute explicitly the inclusive line 
shape of the scalar particle in this toy model, and discuss some features at the next-to-next- 
to- leading order related to gauge invariancc and g^T/M corrections. The generalization 
of the method to final states with detected particles or jets, to non-abclian gauge bosons, 
and to pair production appears not to encounter major new conceptual issues, but requires 
more technical work, which we hope to complete in the future. A brief exposition of the 
effective theory approach discussed in the present paper has been given in [12] . 

The outline of the paper is as follows. In Section 2 we define our toy theory and proceed 
to discuss in a non-technical manner the ingredients of the effective field theory method: 
the momentum scales, the presence of soft and collinear modes, the effective interactions 
and the expression for the forward scattering amphtude in the effective theory. We outline 
a hierarchy of effective theories, in which collinear fields are either kept or integrated out. 
The section concludes with a leading-order calculation of the line shape and a discussion 
of scheme- dependence and matching to the off-resonance cross section. The real virtue 
of the effective theory approach becomes only apparent when one goes beyond leading 
order. A complete next-to-leading order calculation is carried out in Section 3, where we 
compute the two-loop matching of the two-point interaction, and the one-loop matching of 
the production vertex, required at NLO. The NLO result is completed with the calculation 
of the forward scattering amplitude in the effective theory. At the end of Section 3 we 
perform a numerical comparison of the leading-order and next-to-leading-order line shape. 
The formalism is sufficiently general to allow a systematic extension to next-to-next-to- 
leading order. While there is no point of doing this calculation in our toy model, we outline 
the principles of such a calculation in Section 4. Then, as an example of a next-to-next- 
to-leading order contribution, we compute the one-loop short- distance coefficient of the 
F/M-suppressed production-decay operators and illustrate how a gauge-independent result 
is obtained automatically from a combination of self-energy, vertex, and box diagrams. We 
conclude in Section 5. The renormalization constants of the toy theory are collected in an 
Appendix. 
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2 Outline of the method 



In this section we explain the essential features of our approach. We first set up the toy 
field theory and scattering process that we consider in this paper. We then discuss the 
short- and long-distance scales in resonant scattering, the corresponding momentum modes, 
and derive a representation of the scattering amplitude, in which the different scales are 
factorized, allowing for an expansion in F/M. We end this section by giving the two-to-two 
forward scattering amplitude at leading order, and obtain from this the leading-order line 
shape of the resonance. 

2.1 Definition of the model 

We shall consider the resonant production of a massive scalar particle in the scattering of 
two massless Dirac fermions. The scalar and one of the fermions (the "electron" ) are as- 
sumed to be charged under an abelian gauge symmetry, the other fermion (the "neutrino" ) 
is neutral. The charges are supposed to be equal, so as to allow a scalar-clectron-neutrino 
Yukawa coupling. The theory is much simpler than the clcctrowcak Standard Model, but 
we find that all conceptual issues related to the treatment of the scalar resonance can be 
addressed in this toy model. 

We thus assume the Lagrangian 

+ v4>i^x + z/^^^xV' - \{^Hf + ^ct, (2) 

where £ct denotes the counterterm Lagrangian, and the scalar self-interaction is included 
to make the model renormalizable. Fields and parameters are renormalized in the MS 
scheme. The letter M is used to distinguish the MS mass (at scale /i) from the pole mass 
M defined below. Explicit expressions for the counterterm Lagrangian and the rcnormal- 
ization constants are given in the Appendix. The gauge coupling g and Yukawa coupling 
y are assumed to be of the same order. We define ag = g'^/{ATT), ay = (|/|/*)/(47r), and 
use a to refer to them summarily. With no arguments the couplings are evaluated at scale 
yU. In general, we cannot set the scalar self-coupling A to zero, but we would like it to be 
small (for simplicity of the model). Without fine-tuning it is consistent to assume that 
a\/{4:n) = A/(167r^) is of order Q;^/(47r)^, since the leading counterterm is of this order. 
We shall assume this counting in the following. 

We shall consider the totally inclusive cross section in electron-neutrino scattering, 

%) + e-(p)^X, (3) 

as a function of the center-of-mass energy squared s = {q + pY in the vicinity of s ~ M^, 
where we expect an enhancement of the cross section due to the resonant production of 
the charged scalar. By vicinity we mean that s — ~ MF ~ M'^a -C M^. Because the 
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electron is massless, the total cross section is not infrared-safe. The initial-state coUinear 
singularity should be factorized into the electron distribution function. This being under- 
stood, we will usually quote the "partonic" cross section with the singularity regularized 
in dimensional regularization and subtracted minimally. Note that strictly speaking the 
electron distribution function is not defined if the electron is truly massless as assumed in 
(2); however, we may always assume that the electron has a mass much smaller than any 
other scale in the scattering process. We can then neglect the mass in the Lagrangian, 
keeping in mind that it must be reintroduced to regularize the distribution function. 

Our aim is to approximate this totally inclusive line shape of the scalar resonance in 
electron-neutrino scattering systematically in powers of and 

S = — ^ ~ — . (4) 

This cannot be done with standard methods, since there are kinematic enhancements 
proportional to aM^/{s — M^) ~ 1 at every order in the perturbation expansion in a. 
This is the origin of the well-known need for resummation. 



2.2 Effective Lagrangians, effective vertices and representation 
of the scattering amplitude 

In space-time the resonance is produced in the collision with a characteristic formation time 
of order 1/M, lives a much longer time of order l/F, and then decays, again within a short 
time of order 1/M. We therefore expect some kind of factorization between production, 
propagation and decay. This essentially classical space-time picture is corrected, because 
quantum fluctuations exist at all scales. Factorization persists quantum-mechanically, since 
only long-wavelength fluctuations can resolve the details of production and decay (sepa- 
rated by the long time interval l/F) simultaneously. The resulting picture is shown in the 
left-hand graph of Figure 1 below. (The right-hand graph shows the featureless "back- 
ground processes", in which the flnal state of the collision is produced without the reso- 
nance.) Factorization in the presence of quantum corrections can be implemented with the 
effective Lagrangian technique. 



2.2.1 Effective Lagrangian for soft and collinecir interactions 

To obtain the expansion of the line-shape in S, we construct the effective Lagrangian and 
effective vertices for the long-distance contributions to the process. "Hard" effects related 
to quantum fluctuations with momenta /c ~ M are included as coefficient functions by 
matching the effective Lagrangian to the underlying theory. This can be done in ordinary 
weak-coupling perturbation theory, since hard propagators, being off-shell by an amount 
of order M^, do not cause kinematic enhancements proportional to 1/6. Each Feynman 
diagram is broken into a hard and other contributions similar to the non-relativistic ex- 
pansion of Feynman integrals [11]. The hard region corresponds to the Taylor expansion 
of the Feynman integrand in S (counting the loop momentum k ~ M), and goes into 
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the coefficient functions of the effective Lagrangian and effective vertices. The remaining 
regions correspond to momentum configurations, where propagators are near mass-sheU. 
They must be reproduced by the diagrams in the effective theory. 

What are the momentum modes in the effective theory? Parameterize the momentum 
of the scalar particle near resonance by P = Mv + k, where the velocity vector — 1 and 
the residual momentum k scales as M5 ~ F. Then — remains small, of order M5, if 
the scalar interacts with a "soft" fluctuation with momentum of order M5. The effective 
theory therefore contains soft massless modes, and the heavy scalar near mass-shell. (To 
emphasize the fact that the form of the effective theory is very much the same for stable 
and unstable particles, we will also call the heavy scalar near mass-shell a "soft" mode. In 
previous work [10] the term "resonant mode" has been used.) Note that soft interactions 
do not change the velocity vector of the scalar. We therefore define a field with the 

rapid spatial variation e"**^"'^ removed to represent the heavy scalar with momentum near 
Mv. The field (pvix) carries only the residual slow variation in x over distances of order 
1/5. If the heavy particle was a stable heavy quark, the effective Lagrangian that describes 
its interactions with soft quark and gluon fiuctuations would be the heavy quark effective 
Lagrangian [13], which has become a standard tool in heavy quark physics to separate the 
physics on the scales M and Aqcd- The corresponding Lagrangian for an unstable scalar 
is 

£hset = 2M <t>l (iv . - 0, - ^ i^.^^i^r + ^siiPsiJs + X.^ ^Xs + • • • • (5) 

Here Dg = d — igAg denotes the covariant derivative with a soft photon field, and we have 
kept the rclativistic normalization with mass dimension 1 for the non-relativistic scalar field 
(py, hence the factor 2M. Furthermore, ips [xs) stands for the soft electron (neutrino) field. 
Only the leading power in the expansion in 1/M has been given, but there is no difficulty 
in adding further terms. The leading-power Lagrangian contains a single short-distance 
matching coefficient. A, which can be related to the resonance pole position — iMT 
as we discuss below (Section 3.1). In the pole renormalization scheme A equals —iV. 
Since the Lagrangian describes only effects related to the soft scale M5, all dimensionful 
quantities can be assigned a scaling power of 5. (In such assignments, we set M = 1. 
Dimension is restored by inserting the appropriate power of M.) Partial derivatives and 
the soft photon field count as 5. The field (p^ as well as soft massless fermion fields count as 
because their position-space propagator is proportional to ~ 5^. In particular, 
since A is of order F ~ Ma ~ 5, and so is iv ■ Dg, the decay width of the heavy scalar, or, 
more generally, the coefficient A, is a leading-power effect in the Lagrangian, and must be 
included in the scalar propagator. This is the effective theory counterpart of the familiar 
self-energy resummation. 

The effective theory is not yet complete, because the scalar is produced in the scattering 
of energetic {E ~ M) massless particles. To describe these "collinear" modes, we suppose 
that the electron moves with large momentum in the direction of n_, and introduce two 
reference light-like vectors, n±, with — ri^ — and = 2. A collinear momentum 
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is decomposed as 

IX fj, 

= (n+p) ^ + P± + {n-p) ^ , (6) 

where n+p ~ M, n_p ~ M5 and p± ~ M5^/^. The scahng of the small components is 
determined by the interaction of coUinear modes with soft modes, which implies that a 
generic coUinear fluctuation has an off-shellness of order MS. The scahng of the transverse 
component is then fixed by the poles of coUinear propagators. CoUinear modes in an 
effective field theory framework have been discussed previously in B meson decays into 
light energetic mesons [14]. Our set-up is in fact very similar to this case, since the 
production of the scalar resonance represents the "inverse" kinematics of B decay. The 
neutrino plays the role of the weak current, so that we consider the production of a heavy 
particle by scattering an energetic electron on this "current". The effective Lagrangian 
for the interactions of coUinear modes and soft modes is, again at leading power in the 
expansion in S, 

>CscET = 4 (in-D + tip^,- tip^A % - 7 F.^.F^^'' + ..., (7) 

\ m+JJc + 1€ J Z 4 

where V'c denotes the coUinear electron field, which satisfies ^-t/^c — [15, 16]. The 
covariant derivative iDc — id + gAc contains the collincar photon field. The interaction 
with the soft field appears only in in^D = in^d + gn^Ac + gn^Ag. Again there is no 
difficulty in going to higher orders in the expansion; we discuss this below (Section 3.2). 
The 6 power counting is as follows: the components of derivatives on coUinear fields and 
the coUinear photon field scale as the corresponding coUinear momentum components; the 
coUinear fermion field scales as V'c ~ (^^^^ as seen from the propagator. Note that there do 
not exist coUinear interactions with the heavy scalar in the Lagrangian, since the coupling 
of a collincar mode to the on-shcU scalar causes the scalar propagator to become hard. 
These hard effects are already integrated out, and appear as effective production and decay 
vertices (see the following subsection) . Finally, we must add an effective Lagrangian for the 
coUinear neutrino field. Adopting a frame where the electron and neutrino collide head- 
on, the neutrino Lagrangian is obtained from (7) by substituting ipc Xc2, interchanging 
n_ ^ n_|_, and replacing all covariant derivatives by ordinary derivatives. At the level of 
the leading Lagrangian the neutrino is non-interacting. 

To summarize, Feynman diagrams in resonant scattering involve contributions from the 
following three momentum regions: 

hard (h): p ^ M 

soft (s): p~M5 (8) 
coUinear (cl): p±_ ~ MS^I"^, n+p ~ M, n_p ~ M6 

(We note that what we call "soft" here is called "ultrasoft" in [16] and much of the literature 
on soft-coUinear effective theory.) In the general case several types of coUinear modes are 
required, one for each direction defined by energetic particles in the initial and final state. 
For the inclusive line shape we calculate the forward scattering amplitude, so no direction 
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is distinguished in the final state. We then need two sets of coUinear modes, one for the 
direction of the incoming electron, labeled by "cl" (or often simply "c"), the other for 
the direction of the incoming neutrino (labeled "c2"). When the hard fluctuations are 
integrated out, the effective Lagrangian is built from the soft heavy scalar a soft and 
coUinear photon. As and Ad, respectively, a coUinear and a soft electron, ipci and ips, and a 
c2-collinear and soft neutrino, Xc2 and Xs- (In complete generality, we should also introduce 
the coUinear fields ijjc2, ^c2 and Xci, but they appear only in highly suppressed terms, so 
we can ignore them here.) The effective interactions of these modes are described by the 
sum of >Chset and the soft-coUinear effective Lagrangians for the two directions. In order 
to describe resonant scattering with a desired accuracy one has to match the Lagrangian to 
this accuracy, and calculate the scattering amplitude in the effective theory. The relevant 
terms can be identified a priori by applying the S power counting rules. 

2.2.2 Effective vertices 

The effective Lagrangian misses an essential piece of physics. Since £hset does not contain 
coUinear fields, and since £scet does not contain the heavy scalar field, the two can interact 
only through soft modes. In particular, there is no vertex in the Lagrangian that allows 
the production of the scalar in the scattering of energetic particles. 

Such vertices cannot be included in the effective Lagrangian for soft and coUinear 
terms as interaction terms, because they contribute to the scattering matrix element in the 
effective theory only in a very specific pattern. For instance, a vertex with fields 0i,V'ciXc2 
can occur exactly twice, once for the production of the scalar, and once for its decay. 
Multiple insertions are not compatible with the kinematic restrictions on the process, which 
allow only one nearly on-shell scalar line in any diagram by energy conservation. Note that 
a diagram in the full theory may of course have many internal scalar lines. However, these 
will generally be far off-shell ("hard") and these hard effects are included in the coefficient 
functions of the Lagrangian and effective vertices. What we are discussing here is the 
structure of possible diagrams/scattering processes after hard effects are integrated out. 

In more technical terms the interaction of two coUinear modes with opposite directions, 
here the electron and the neutrino, produces a hard momentum configuration with invariant 
mass squared of order M^, and therefore cannot be included in the effective Lagrangian. 
What we actually need is that the momenta of the oppositely moving coUinear modes are 
pre-arranged (by the experimenter who sets up the beam energy to be near resonance) to 
produce a configuration with invariant mass squared equal to within a small amount 
of order M'^d ~ s — M^. To account for these configurations we introduce production 
and decay vertices ("effective vertices"). The leading-order vertex is simply the original 
Yukawa coupling expressed in terms of the fields in the effective theory, 

J(x)=e-^^''-^?/[(/)„VSciXc2](x). (9) 

In higher order in a or 5 a larger set of operators - not necessarily local - is generated by 
integrating out hard quantum fluctuations and the coefficient function will receive correc- 
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Figure 1: Reduced diagram topologies in 2 — > 2 scattering near resonance. Left: 
resonant scattering. Right: non-resonant scattering. See text for explanations. 

tions. (This will be seen in more detail in Section 3.3.) However, these operators have in 
common that they contain the field (f)^ exactly once. 

The scattering may also occur without the production of the scalar near its mass-shell. 
In our toy theory this still requires an intermediate scalar line, since the neutrino has only 
Yukawa interactions. The scalar may be off-shell, because the electron has radiated an 
energetic (hard or collinear) photon before it hits the neutrino. In this case the invariant 
mass of the colliding electron-neutrino system is of order but not near M^, producing 
a non-resonant scalar. In the effective theory this process is represented by production- 
decay operators, which do not contain (p^ fields. The simplest operator of this type has the 
coupling and field structure 

T{x) ~ \y\^ [{^ciXc2){Xc2M]{x). (10) 

In general, non-resonant scattering includes all "background processes" , which produce one 
of the final states under consideration. 

With these preparations we obtain the following formula for the forward 2 — > 2 scat- 
tering amplitude in the effective theory: 




where in the first term "T" denotes time-ordering, and where the sums extend over the 
sets of effective vertices of the two types, J or T. The matrix elements on the right-hand 
side are evaluated with the effective Lagrangian for soft and collinear fields. This result is 
depicted in graphical form in Figure 1. The "resonant" scattering term on the left-hand 
side involves a production and decay vertex, and collinear interactions at the electron 
legs, which connect to the production and decay vertex. Soft (and only these) fields can 
connect to all parts of the diagram. The double line denotes the heavy resonance, solid lines 
collinear modes and dashed fines soft modes. A single sohd (dashed) line may represent 
a collection of several collinear (soft) fields, depending on the structure of the effective 
interaction. Collinear fields cannot be exchanged from the left to the right of a resonant 
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diagram, since this would cause the scalar to be off-shell. The corresponding configurations 
are included in the "non-resonant" scattering topology on the right-hand side of the figure. 
This topology does not involve a resonant heavy scalar, and both soft and coUinear fields 
can be exchanged across the diagram. 

The various terms in (11) can be ordered according to their scaling with 5 and a, thus 
allowing for a determination of the terms needed for a calculation to some accuracy before 
performing any computation. We will now estimate the importance of the resonant and 
non- resonant contributions by a simple power counting argument, taking as representatives 
the leading interactions of each type. Combining the scaling of the fields we find J ~ 
from (9). Furthermore x is soft, since x represents the distance of order 1/S over 
which the resonant scalar propagates, hence d'^x in (11) counts 1/5^. The normalization of 
states implies that an external coUinear particle counts as 1/5^^^, so |Pe) ~ 1/S. Putting 
everything together, we obtain 

a^ilmT^^^i|L, (12) 

which is the expected result. The leading non- resonant scattering operator T in (10) scales 
as The matrix element (Pe|iT(0)|Pe) is non- vanishing only at order g'^, since the 

overlap of T with the external states requires the emission of an energetic photon. Adding 
the counting of the states, we find that the contribution to T is of order g'^\y^\. This 
is suppressed by a factor g'^V/M relative to (12), hence establishing that non-resonant 
topologies contribute only at next-to-next-to-leading order in the combined expansion in 
a and 5. We shall discuss these effects in some detail in Section 4.2. 

In this discussion of the non-rcsonant contribution wc assumed that the {ipciXc2) op- 
erator in T describes an electron-neutrino configuration with invariant mass far from M^. 
We shall see below that the applications of the equation of motion to bring the scalar ef- 
fective Lagrangian into its canonical form generate a local four-fermion operator identical 
to T, but with (■0ciXc2) describing an electron-neutrino configuration with invariant mass 
near M^. This arises because the resonant propagator is canceled by the rearrangements 
implied by the equation of motion. In this case T has a non-vanishing tree-level matrix 
element and contributes already at next-to-leading order. 



2.2.3 Soft and coUinear fluctuations of coUinear modes 

The complication mentioned in the previous paragraph, as well as the complication of 
excluding "by hand" coUinear loop contributions to the matrix element of T | J^(0) J(x)|, 
which connect the initial and final state, can be avoided, if we distinguish modes that differ 
from the external coUinear momenta only by a soft momentum from those which differ by 
a collincar momentum. Both modes have virtuality Md^^"^ and have been called "coUinear" 
up to now. However, the former can be assigned momentum Mn_/2 -|- k, where k is soft 
and the large coUinear component Mn_/2 is fixed. Only the latter represent coUinear 
fiuctuations associated with coUinear loop momenta. 
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Just as for the heavy scalar with momentum close to Mv, it is useful to extract the 
fixed large momentum from the field, and to define 



(13) 



for soft fluctuations around the external momentum. {V+ projects on the positive frequency 
part of ipci-) The new field ■(/;„_ (x) has only variations over distances of order 1/6. For the 
collinear fiuctuations we keep the field ipci {x) . Similarly we define the field Xn+ {x) . 

The collinear effective Lagrangian now assumes a more complicated expression, since we 
distinguish the interactions of the two types of collinear fields. However, it is now possible 
to add the effective vertices as interaction terms in the Lagrangian, since the kinematic 
distinction of configurations with virtualities of order and near is implemented at the 
level of the fields. In particular, the effective vertex operator J{x) is now given by 



but there is no interaction vertex of this form with ijjn_ replaced by i/jd- For the four- 
fermion operators, we must distinguish {'il'ciXn+){Xn+'4'ci) from {'il'n-Xn+)iXn+'4'n-) ■ Since 
the external coUinear state can be created or destroyed only by the softly fluctuating 
collinear fields and their conjugates, the latter operator has a tree-level matrix element, 
but the former does not, because it requires interaction vertices of the type ■ipn-i'ci^ci from 
the collinear Lagrangian to generate a non-vanishing overlap with the external state. 

2.2.4 Integrating out collinear fluctuations 

The formalism developed up to this point contains propagating soft and collinear modes 
in the effective theory. We can go one step further and integrate out collinear fluctuations, 
that is we integrate out loops with momenta that scale as collinear momenta. This can 
be done perturbatively as long as the couplings at the scale M5^/^ are small. We now 
sketch the structure of the resultant effective theory. At the level of the calculations in the 
remainder of the paper there is in fact no difference between the two descriptions before 
and after integrating out collinear loops, since the collinear loop integrals we encounter all 
vanish in dimensional regularization. 

Technically we introduce two distinct collinear fields as described above and eliminate 
the fields that describe collinear fluctuations from the effective Lagrangian. With only soft 
modes remaining in the effective theory, the effective Lagrangian is constructed from the 
soft electron and photon held, and the flelds 0^, i^n--, and Xn+i all scahng as 5^/^. The 
scattering amplitude can be computed from 



where £hset describes the soft interactions of the heavy scalar as before, and 

are bilinear in the flelds ipn^ix) and Xn+{x), respectively. The soft-collinear Lagrangian 

>CscET is non-local due to the presence of l/(in_|_9), but these non-localities disappear. 



y[(t>v'^n-Xn+]{x), 



(14) 



^ — -^HSET + ^- + ^+ + -^int) 



(15) 
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when coUinear fluctuations are integrated out, and the large component of the energetic 
fields is "frozen" to M/2 (modulo soft variations). For instance, consider the leading-power 
SCET Lagrangian (7) for the electron field, which matches at tree level to the expression 

^n_(^n_I^. + J^)^^n_ (16) 

in the new effective theory. The term with the transverse derivative squared is already 
suppressed by one factor of 6 compared to the leading term, since derivatives on ipn^ix) 
count as 5. When coUinear quantum fluctuations are included, the set of allowed operators 
consists of operators non-local in l/(m_9), since the n_ component of the external soft 
momenta of a coUinear loop are of the same order as the n_ component of the coUinear 
momentum. These new non-localities, appearing at higher orders in a, correspond to the 
convolution of a "jet factor" with the remaining soft matrix element. Similar considerations 
apply to the neutrino field. The fields for the soft fluctuations of the energetic particles 
now behave in many ways similar to the field for a static heavy particle. In both cases 
the virtuality is M'^S, and the interaction with soft modes (of virtuahty M'^S'^) kicks the 
particle off its classical trajectory only by a negligible amount; hence the fixed vector that 
labels the field. (The discussion here is reminiscent of "large-energy effective theory" [17], 
except that we have in mind a situation where hard and coUinear modes are eliminated.) 

As mentioned above, the production and decay vertices can be implemented as interac- 
tion terms Cmt in the Lagrangian, because the kinematic restriction that prevented us from 
adding generic interactions of coUinear modes with different directions to the Lagrangian 
is now implemented at the level of the fields: by construction the scattering of the elec- 
tron described by (x) and the neutrino described by Xn+{x) produces a mode with 
momentum Mv + k with k soft. The most general terms allowed for a 2 — > 2 scattering 
process are of the form (f)v'4^n-Xn+ {or its hermitian conjugate) or (^n_Xn+) {Xn+'4'n-) plus 
additional soft fields. Since none of the fields 0^, and Xn+ is created or destroyed by 
the other terms in the Lagrangian, the terms in Ci^t can only contribute in a topology that 
corresponds to a classical scattering process of the energetic particles and the heavy scalar. 
(In particular closed loops of only 0^, ipn. and Xn+ lines vanish due to the structure of 
the poles of the propagators.) These are the two topologies shown in Figure 1 but with 
the coUinear subgraphs eliminated and interpreted as parts of the production and decay 
vertices. The power counting goes as follows: since / d'^X(pvipn_Xn+ ^ S^^'^, each insertion 
of this interaction vertex gives a factor of 5^/2^ (E ven though the effective theory contains 
only a single scale, we must count powers of 6 rather than dimensions due to our uncon- 
ventional 2M normalization of the non-relativistic scalar kinetic term.) The insertion of 
a four-fermion operator carries a 5^ suppression factor, which results in a 5 suppression 
relative to the resonant contribution, not counting factors of coupling constants. 

2.3 Leading- order line shape 

We now perform the (trivial) leading-order calculation of the inclusive line shape of the 
resonance. From the foregoing it follows that we need only the leading-order Lagrangian (5) 
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together with the leading vertex J in (9). The calculation consists of two steps: matching 
the effective Lagrangian to the "full" theory to leading order, and the calculation of the 
forward scattering amplitude in the effective theory. At the end, we use the optical theorem 
to obtain the line shape. 



2.3.1 Matching 

Since ivd is of order 5 on soft fields, we see from (5) that we need A to order 5 ~ a. 
This implies that A must be chosen so as to reproduce the (renormalized) full theory 
propagator i/{p'^ — — 11 (p^)) near the resonance pole with this accuracy. We discuss 
the matching condition in more technical terms in Section 3.1. At leading order, one finds 
A*^^) = n^^^(M^)/M, where n'^^^(M^) is the renormalized one-loop self-energy evaluated at 
= M\ 

All calculations will be done with dimensional regularization with ci = 4 — 2e, where 
the loop integration measure is defined as 

^ ^ V 4vr ; {2'kY ^ ' 

With this convention the MS scheme corresponds to subtracting all poles in e minimally. 
In our toy model a straightforward computation of four diagrams 




n« = -----W----- + ---/ \— + ---) V— + X (18) 



gives 

A^ ay f , \ a, f , \ 

^ = -^ 21n— -4-2m -31n— + 7 19 

M 47l\ fl^ J 47l\ /i^ / 

in the MS scheme. (Note that the one-loop diagram proportional to the scalar self-coupling 
A counts as and therefore does not contribute to 11*^^).) 

The renormalization scheme dependence of A is related to the mass renormalization 
convention in the underlying theory, if we choose the parameter M in the effective theory 
to be the renormalized mass of the full theory. The effective theory formalism does not 
rely on a specific renormalization convention. We will use the MS convention and the pole 
mass convention to illustrate this point. The pole mass M is defined through the location 
of the singularity, 

s = M^- iMT (20) 

in the scalar propagator. Writing M = M + 6M, and using A^^^ = n(^)(M^)/M, the 
relations 

5M« = ^ReA«, r« = -ImA« (21) 

hold at the one-loop order. Note that M = M{fi), but we do not indicate the scale- 
dependence explicitly. 
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2.3.2 Forward scattering amplitude 

The forward scattering amplitude is given at leading order in the effective theory (not 
distinguishing explicitly between coUinear and soft fluctuations of the collinear fields) by 

= i^yy* \u{v)v{q)] — — \v{q)u{v)\, (22) 

where the double line denotes the resonant propagator defined by the effective Lagrangian 
-Chset, and k = p + q — Mv. We take v = (1, 0) and the external momenta p = 
v^/2 (1,0, 0,-1) = ^/s/2n-, q = v^/2 (1, 0, 0, 1) = v^/2r2+, so v-k = ^s-M. Note that 
by construction the effective propagator includes the geometric series of all one-loop self- 
energy insertions, evaluated at threshold according to the definition of A*^^). Performing 
the polarization average the forward scattering amplitude reads 

= yyl^ f23) 

4M(v/^-M- A{i)/2)' ^ ' 

Note that because we used spinors u{p) etc. for the external states, T^^^^ contains next- 
to-leading order terms from the factor s = + — M'^] in the numerator. For a strict 
expansion it would be more appropriate to use the spinors M(Mn_/2) etc., in which case 
T*^°) would be strictly leading order, but in practice it may be convenient not to perform 
this trivial expansion in order to reduce the number of terms. 

2.3.3 Inclusive line shape 

The inclusive line shape is related to the forward scattering amplitude by the optical 
theorem, which gives 



s 4M (v^i- [M + 5M«])2 + r«74 

= . J (24) 

4 {^-[M + 6Mm]y + {ayM/4y 

where we used (21). The line shape has the form of a Breit-Wigner distribution in the 
center-of-mass energy -^i. This is in fact the universal line shape that appears in the 
leading-order approximation to any resonant scattering process in our toy model, but 
the line shapes will depend on the final state beyond this approximation. (It is a matter of 
convention whether we divide by s to obtain the cross section from the forward scattering 
amplitude, or expand this factor around M^. Here we keep 1/s unexpanded, so that the 
factor of s in the numerator of (23) is canceled, and the leading-order line shape assumes 
an exact Breit-Wigner distribution.) 
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In the presence of fields for unstable particles unitarity and the optical theorem apply 
to the scattering matrix defined on the Hilbert space of stable particle in- and out-states 
only, since the resonance does not correspond to an asymptotic particle state [18]. The 
optical theorem can be interpreted as taking the sum over all "cuts" of a diagram, where 
now the line for the unstable particle propagator must not be cut. For the leading order 
diagram in (22) this means that cutting the effective propagator represents the sum over all 
possible cuts of the one-loop self-energy insertions implicitly contained in the double-line 
propagator. This corresponds to cutting the second diagram in (18), which is responsible 
for the leading-order decay width T^^^ of the scalar. 

We remark that the effective Lagrangian is not real, since the hard coefficient functions 
have imaginary parts, but this does not lead to a non-unitary time evolution. Consider, 
for instance, the diagram in (22). Unitarity requires an amplitude that corresponds to the 
"square" of this diagram, but the effective theory does not contain diagrams with closed 
electron-neutrino loops. Rather the corresponding (short-distance) effect is included in 
the complex coefficient function. Matching the effective theory to the underlying theory 
automatically guarantees that the combination of effective theory diagrams and complex 
couplings reproduces the unitary time-evolution to the desired accuracy. 

The calculation shows that if the result is represented in two different mass renormal- 
ization schemes, the masses must be related to one-loop accuracy for a consistent change of 
conventions. For instance, if the pole mass M — 100 GeV is known from elsewhere, the MS 
mass ioT fjL — M is M — M — SM'^^^ — 98.8 GeV, where we use ay — Ug — 0.1 to illustrate 
the numerics. From (24) we see that the denominator contains ^/s — M in any scheme to 
one-loop accuracy, but there is a residual scheme dependence (and scale dependence, since 
M depends on /i) due to the width ayM/2. This is a NLO effect, which will be reduced 
by a higher-order calculation, as will be the dependence on the renormalization convention 
of the coupling constants. The line shape in the pole and MS scheme is shown in the left 
panel of Figure 2. 

The effective theory calculation (24) is valid only near resonance when ^Js — M is small, 
since the expansion inherent to the construction neglects corrections of this order. In other 
words, near resonance the corrections are of order a or 5, but the second quantity becomes 
large far off the resonance. To obtain an adequate description of the line shape in the 
entire kinematical range, the effective theory calculation must be matched to a full theory 
calculation off resonance. The point to note here is that the full theory calculation of 
the forward scattering amplitude can be done with ordinary perturbation theory, since no 
kinematic enhancements invalidate it off resonance. In particular, the optical theorem is 
applied to intermediate states including the scalar particle as if it were stable. 

The off-resonance cross section is only of order compared to order one near resonance. 
The result takes the form 

with 

fy{z) = 1, 
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Figure 2: Left panel: line shape (in GeV^^) in the pole (solid) and MS scheme 
(dashed) as a function of the center-of-mass energy (in GeV). Right panel: Leading- 
order line shape (in GeV~^) as a function of center-of-mass energy (in GcV) in the 
effective theory (solid) and the cross section off resonance in the full theory (dashed). 
The thick grey curve shows the leading-order line-shape with the two curves matched 
appropriately. 



m = (l-^)|(l + ^^)ln^LiL_2^|. (26) 

The "partonic" cross section has a coUinear divergence, which has been subtracted accord- 
ing to the MS prescription. The physical cross section follows from convoluting the above 
result with an electron distribution function (in the electron) that depends on the subtrac- 
tion scale u. The argument of the logarithm suggests that the collinear factorization scale 
should be taken to be of order (s — M^)/M ~ MS. In the right panel of Figure 2 we show 
the leading-order computation in the effective theory, and in the full theory off resonance, 
evaluated at the running scale u — ^/s (1 — z). The effective theory calculation correctly 
represents the line shape near the resonance with relative accuracy a. The peak height is 
of order 1. As the calculation is extrapolated off resonance the cross section decreases and 
becomes of order a^. However, in this region the relative error on the effective theory cal- 
culation is of order unity. The full theory calculation correctly represents the off-resonance 
cross section of order with relative accuracy a. As the calculation is extrapolated to- 
wards the resonance, the relative error becomes of order unity. (In fact, the cross section 
diverges at s = and the error becomes arbitrarily large.) The two calculations must 
be matched in the intermediate region where 5 ~ (y^ — M)/M is small enough for both 
calculations to be valid. The existence of this intermediate region is indicated in the figure 
by the vertical lines to the left and the right of the resonance peak. We also verify from (24) 
and (25) that the two calculations agree analytically in this intermediate region. (For com- 
pleteness we mention that the partonic cross section proportional to ayS{s — M^) should 
be included in the convolution with the electron distribution function, which ensures that 
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the off-resonance cross section is i/-independent, and sums large logarithms of ^/s over the 
electron mass. Since this procedure is standard, we presented the partonic cross section at 
the running scale, at which no large logarithms occur.) 

3 The line shape at next-to- leading order 

Following the formalism set up in the previous section we now discuss the matching calcu- 
lations in more detail. We construct the effective Lagrangian by matching on-shell Green 
functions and present the explicit form of all terms of the effective Lagrangian that are 
needed for the evaluation of the line shape at next-to-leading order. At the end of this 
section we also compute the scattering matrix element in the effective theory to next-to- 
leading order and obtain the complete NLO line shape. 

The matching procedure involves three steps. Compute the renormalized on-shell Green 
functions in the full theory up to the required order in a. Evaluate the same quantity in 
the effective theory. Determine the hard matching coefficient so that the two calculations 
agree within a specified accuracy. The Green functions in the full theory are calculated 
using conventional weak-coupling perturbation theory. In general, we need an infrared reg- 
ularization, and the same regularization must be used when computing the corresponding 
quantity in the effective theory. The result for the matching coefficient is independent of 
the IR regularization. The calculations become particularly simple with dimensional reg- 
ularization, since the effective theory loop diagrams vanish in this scheme. This is simply 
a consequence of the fact that the matching is done on-shell, which renders the corre- 
sponding integrals scaleless. However, the matching coefficients depend on a factorization 
scheme, since the Green functions in the effective theory also exhibit ultraviolet singular- 
ities. Thus they have to be UV-subtracted and the choice of the renormalization scheme 
for the effective theory amounts to the choice of the factorization scheme. 

In practice, we simplify the matching procedure by computing directly the hard part 
of the renormalized on-shell Green functions in the full theory, using the method of re- 
gions [11]. This also applies to the Z-factors that multiply the amputated on-shell Green 
function. Even though the method of regions forces us to use dimensional regularization 
to define the factorization scheme, we are still free to use another infrared regulator, such 
as a photon mass, in which case the loop diagrams in the effective theory are generally 
non- vanishing. However, the result of the hard part is independent of the infrared regu- 
larization of the full theory, but it contains infrared 1/e poles due to the separation into 
hard and soft contributions. These singularities coincide with the above-mentioned UV 
singularities in the effective theory We choose to subtract them minimally and thereby 
define the factorization scheme and the renormalization scheme for the effective theory. 

3.1 Matching /^hset 

For the calculation of the line shape at leading order it is sufficient to consider the first 
term of £hset [Eq. (5)], with the matching coefficient A known to order a ~ 5. At next-to- 
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leading order we will need all terms of order {a^, a5, 5^} x 01^0^,. Furthermore, at NLO there 
will also be soft photon loop diagrams in the effective theory, so we will need to include 

the corresponding terms in the Lagrangian. Since these diagrams will be suppressed by a 
factor a with respect to the leading-order diagram, the soft photon vertices are only needed 
at leading order. 



3.1.1 The bilinecir terms for the unstable field 

The bilinear terms of the effective Lagrangian are determined by the dispersion relation 
for the free-particle states. The field describes a heavy particle near mass-shell, so we 
need the equivalent of the heavy-quark effective Lagrangian for a scalar field. The only 
complication is that the scalar is unstable, so the notion of a free-particle state is not really 
defined. In this case the bilinear terms are constructed so as to reproduce the two-point 
function near the resonance pole. 

In the underlying theory the full renormalized propagator for the unstable particle is 
given by i(s — — n(s))~^, where — in(s) corresponds to the amputated IPI graphs, 
including counterterms. Recalling P — Mv + k the propagator near the resonance pole 
can be written as 

P'^-s 2Mvk + k"^ - (s - M^y 

where we denoted the complex pole of the propagator by s and the residue at the pole by 
R(j). Defining 

A.i^, (28) 

we solve = s in the form 

A A^ -I- Ak^ 

vk^-M + JM'^ + MA-kl^- -^ + 0(6'^), (29) 

V ' 2 8M 

where for any vector ay = a'* — {va) v^. The second solution has vk M and is irrelevant. 
Prom this we obtain the effective Lagrangian 

£^^ = 2M(/)t [iv.Ds-^^(t>v + 2M(t>l + + . (30) 

The terms in the second bracket are suppressed by one factor of 5 relative to the leading 
terms. Several comments are in order: 

(i) The field 0^ is a pure destruction field. Even if the original scalar field had been 
neutral (real), the corresponding effective scalar field would have been complex, because 
it contains only the destruction part of the relativistic field. We keep a factor 2M in the 
kinetic terms to preserve the canonical mass dimension 1 for a scalar field. 

(ii) Since the full Lagrangian is gauge invariant and so is the separation into hard 
and soft contributions, the effective Lagrangian is also gauge invariant. Therefore, the 
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interactions of the unstable field with soft photons can be taken into account by replacing 
ordinary derivatives by covariant derivatives. This has been done in (30). 

(iii) The quantity A plays the role of a matching coefficient. For a stable particle A 
vanishes in the pole mass scheme and is referred to as "residual mass" in a general scheme. 
In the unstable case A is complex and always non- vanishing. It is A that prevents the 
propagators of the effective theory diagrams from becoming singular. For the interpre- 
tation of A as a matching coefficient to be consistent, A must be given entirely by hard 
fluctuations. This is related to the fact that the location of the pole of the propagator (27) 
is infrared-flnite. Alternatively, we may note that quantum corrections to the effective the- 
ory propagator do not modify the location of the pole to any order in perturbation theory, 
so the full self-mass must already be contained in the coefficient function A. Further, we 
note that A is gauge-independent, because s and M are gauge- independent. 

(iv) Computing the propagator to all orders in S with the effective Lagrangian (30) 
does not reproduce (27). Instead we obtain 

2M (ivP) - ^Jm'^ + MA - P2 _ s 

near the resonance pole, where we used that the coefficients of the higher-order kinetic 
terms in the effective Lagrangian renormalizc in the same way as the leading term due to 
Lorentz invariance. Besides the different standard residue factors, there also appears the 
factor 

^-^(l + ^^^^y'^l + 0(a,5), (32) 

which differs from 1 even at tree level (A = 0). 

This is due to the fact that (up to the factor 2M) we use the standard form of a 
non-relativistic Lagrangian, which does not reproduce the normalization of a relativistic 
field. This is not a problem, because the different normalization is taken into account in 
matching calculations. Whenever we compute an amputated Green function in the effective 
theory we multiply every external (py-lvae by the additional wave-function normalization 
factor ti7~^/2. This will be important in Section 3.3, where we discuss the matching of the 
production (decay) vertex. 

3.1.2 Two-loop computation of A 

The matching coefficient A is of order M5 with 6 defined in (4) and will be computed from 
the perturbative expansion of the hard part of the self-energy only. We write the expansion 
in the form 

n^(s) ^M^^^'n^'^''), (33) 

k,l 

where it is understood that n^'^'^) ~ a^. Since the unstable field couples to massless particles 
the full self-energy is not analytic at M^. In particular, for our model dn(s)/ds|^^jg2 has 
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an infrared singularity. However, only the hard part of the self-energy enters (33). This 
part is constructed from the Taylor expansion of the Feynman integrand in 5 and hence 
analytic at M^, so the Taylor expansion (33) is also well-defined. 

The position of the pole of the propagator, s, and the hard part of the residue at the 
pole, Rh<i>, can be expressed in terms of the n^'^''^ Up to the third order we find 



^ ^ i + n(i,o)^n(2,o)^n(i,i)n(i,o) (34) 

+ n(^'°) + n^^'^^n^^'") + n(i'i)n(^'°) + [n^^'^^j^n^^-") -1- n^^-^^p^^-")]^ + .... 
Rj,l = 1 - n(i.i) - n(2.i) - 2n(i'2)n(^'°) (35) 
_ n(3.i) _ 2n(2'2)n(^'°) - 2n(^'2)n(^'°) - 2n(^'2)n(^'^)n(^'°) - 3n(^'^)[n(^'°)]2 + .... 

The first equation allows us to express the matching coefficient A [Eq. (28)] in terms of 
the hard part of the self-energy. Note that even though only the hard part of 11 appears in 

(34) , s is the physical pole location. In other words, there is no difference whether we use 
the full self-energy or only its hard part in (34), since A receives no soft contributions as 
discussed above. The situation is different for the residue, since its hard part, as defined in 

(35) , does not in general coincide with the residue of the full propagator R^ that appears 
in the LSZ reduction formula and in (27). i?^ is given by an equation similar to (35) where 
the hard part is replaced by the full self-energy The relation reads R^^ = Rh<j>RcSi)>- The 
residues i?<^ and -Rcflf^ in the full and effective theory arc IR divergent and depend on the IR 
regulator, but the hard part of the residue is independent of this regularization. However, 
if dimensional regularization is adopted for the IR singularities, ReSip = 1 (since all loops 
in the effective theory vanish), and the full and hard residues coincide. 

We now turn to the explicit calculation of the matching coefficient 

A = J2 A^'^ = M n(^'°) + M (n(2.o) + n(i'i)n(i'°)) + . . . (36) 



with A^'^ = O (a*). We need H(^'°^ for the calculation of the matching coefficient at leading 
order, while we need H^^'^^ and H^^'^^ at NLO. Defining 

"^^4^' "^"4^' "^^16^' (^^^ 



we find 



n(i'O) = s(^) (^-^-4-8e+^-2i7r-4i7re^ (3J 

+ %(f)~7- + r + i5. + 4^V.£)-e 

na,i) = ^ l---2 + ^-4e-2t7T-2t7:e]-6^;\ (39) 
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where and 5^^^ refer to the order a part of the counterterms given in the Appendix. 
The order e terms will be needed only in (41) below. 



The two-loop self-energy has contributions proportional to ctyCig (4 diagrams) 



a 



y 



diagrams) and (10 diagrams), not counting counterterm diagrams, as well as a contri- 
bution proportional to ax, which by assumption is of the same order. The result reads 
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with the counterterms given by 



n, 



(2) 
ct 



+ 



9M2 



In (41) we indicated the power of a in the counterterms and denoted by 11^^''') and Il^g'^^ the 
contribution to the one-loop self-energy proportional to ay and ag respectively, not includ- 
ing counterterms. We adopt the MS renormalization scheme and refer to the Appendix for 
the explicit expressions of the counterterms appearing in (41). The two-loop integrals have 
been computed with standard techniques. We checked some of our results with recurrence 
relations from [19], and two- loop master integrals from [20]. Massive one-loop integrals 
(needed here and later) were in part computed with the method described in [21]. 
The matching coefficients A*^^) and A^^^ defined in (36) read in the MS-scheme 
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We have computed these coefficients in an arbitrary covariant gauge and checked their 
gauge-parameter independence. Note that 11*^^'^) and 11*^^'°^ in (39) and (40) have 1/e 
infrared poles after MS renormahzation, which must be associated with the factorization 
of the hard and soft contributions. For s and, therefore, for the matching coefficient A 
these singularities cancel as they should. On the other hand, the residue at the pole, R^p, 
has an infrared singularity in the MS-scheme. (In the on-shell scheme i?^ = 1 by definition, 
but in this scheme the field renormahzation factor exhibits the infrared singularity.) 

3.1.3 Including soft photons and fermions 

To complete the soft Lagrangian £hset for the NLO computation the bilinear terms of the 
unstable field, £<^0, have to be supplemented by the kinetic terms of the soft photon and 
fermions. Due to gauge invariance, these terms will also contain the interaction terms of a 
soft photon with a soft charged fermion. 

The leading operators scale with 5^ and are given by the gauge-invariant kinetic terms 
ipsi pgi/jg, XsiPXs and —1/4 Fs^uFjf'" ■ By convention the kinetic terms are canonically 
normalized, and their coefficients do not receive corrections from hard loops. The effective 
Lagrangian contains an infinite number of higher-dimension operators such as ipsi^sXsXs 
(which scales with S^). At NLO we need the operators that scale with S^. However, since 
we assume weak coupling with g ~ 5^^^, at NLO we only need the operators that have 
non-zero coefficients as g 0. There are no such operators, so we conclude that at NLO 
the only relevant soft photon-fermion interaction is contained in i/jgips'tps- No further 
calculation is required to match £hset at NLO. 

Putting everything together we obtain for £hset including all terms needed for the 
NLO line shape in the weak-coupling limit 

£hset = '2M(Pl (^iv ■ D, - ^ j (f>, + 2M(f>l 

- ^Fs^^F^ + i^stPsi^s + Xst^Xs (44) 

with A^-*^) and A^^^ given in (42) and (43) respectively. The first term of (44) gives rise 
to the unstable particle propagator used in the calculation of the leading-order amplitude 
(22). All other terms are needed at next-to-leading order only. Note that in writing the 
Lagrangian (44) we have expanded the matching coefficient A in powers of a, so that 
the unstable particle propagator strictly counts as 1/5 ~ 1/a. Alternatively, we could sum 
higher-order corrections to A back into the propagator, although this would not correspond 
to a strict expansion of the scattering amplitude in 5 and a. We will investigate the effect 
of including additional terms in the propagator in the numerical analysis below. 

3.2 Matching jCscet 

The leading Lagrangian for the collinear fields has been given in (7). In the position space 
formulation of soft-coUinear effective theory the Lagrangian has been worked out to order 
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S in [16]. (We note again that what is called "soft" here has been called "ultrasoft" in 
[16].) When expanded systematically in powers of S^^'^, the interaction vertices are rather 
complicated. Our task of writing down all the terms needed for the NLO computation 
of the line shape is again simplified by the fact that we adopt the weak-coupling limit 
g ~ 5^/^. Inspection of the expressions in [16] shows that once again there are no relevant 
operators at order 5^/^ and that survive as ^ 0. All the vertices needed for coUinear 
loop corrections to the LO line shape are contained in the leading-order Lagrangian 

V^c (in-D + ilp^^- ilp^A % Vc- (45) 

The Feynman rules for the three-point vertices 'ipcipc^c and ipcipc^s can be read off from 
this expression. We may note that the inverse covariant derivative can be written in terms 
of Wilson lines 

{in+Dc + ie)-^ = Wc {in+d + ie)'^ (46) 

where Wc is 

Wc{x) = exp J dsn+Ac{x + sn+)^ . (47) 

The term (45) therefore contains vertices with any number of Ac fields, which are all 
leading power in 5, but suppressed by gauge coupling factors g, so they are not needed for 
the computations of the NLO line shape. 

Interactions with soft photons are contained in 'ipciri-D ipc- In [16] coUinear-soft in- 
teraction vertices are light-cone multipole-expanded, and the covariant derivative reduces 
to in_D = in_d + gn_Ac{x) + gn_As{x_) with ,t_ = (n^x) n_/2 at leading power. This 
should be done here only for the coupling of soft fields to the fields describing coUinear 
fluctuations, but not for the coupling to coUinear flelds with soft fluctuations, since only 
for the former fields is the transverse momentum large compared to the transverse mo- 
mentum of soft fields. This affects the way momentum conservation is implemented at the 
corresponding interaction vertices. 

3.3 Production (decay) vertex 

The computation of the LO forward scattering amplitude (line shape) requires only the 
vertex (/)v4'ciXc2 + h.c. at tree level. At NLO we therefore need the one-loop coefficient 
function of this vertex, and production or/and decay vertices suppressed by one factor of 
5 (without loop corrections). 

3.3.1 Power-suppressed vertices 

The power-suppressed vertices with field content 4>vi^ciXc2 can be obtained to order 5 from 
the result for the heavy-to-light decay current in [16]. However, the result relevant for 
the forward scattering amplitude can be obtained in a simpler way, since we only wish to 
reproduce the tree 2^2 scattering diagram (22) to better accuracy. 
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For tree-level matching we can set n(s) = 0. Leaving out the external spinors and 
setting P — Mv + k, we obtain by expansion in the small momentum k 



iy —iy* = iy — x iy* + iy — ik^ — iy* + i-^ — h 0(S). (48) 

The first term on the right-hand side is accounted for by the leading-order term in Z^hset 
together with the leading-order production (decay) vertex. The second term involves one 
insertion of the kinetic energy correction (^-D^y) in (44). (For head-on collisions this term 
vanishes at tree level, because fcj = 0.) The third term is not reproduced by any term in 
the Lagrangians we discussed so far. Since the intermediate unstable scalar propagator has 
been canceled by the expansion, we associate this term with the production-decay operator 

T=fC-ii^ciXc2){Xc2M. (49) 
4M2 

It should be noted that this operator contributes at tree level contrary to the off-shell 
scattering topologies in the full theory discussed in the previous section. When we dis- 
tinguish coUinear fields with coUinear fiuctuations from those with soft fiuctuations, the 
operator discussed here appears with field content (■0n_Xn+)(Xn+'*/'n_); since each fermion 
pair describes a configuration with invariant mass close to M^. This must be distinguished 
from the off-shell configurations produced in the scattering of a coUinear electron fluctu- 
ation off the neutrino. This scattering process is represented at tree-level by a non-local 
four-fermion operator, but it does not contribute at tree-level to the forward scattering 
amplitude in the effective theory, resulting in a NNLO contribution. 



3.3.2 Hcird loop correction 

We now turn to the one-loop matching of the coefficient function C of the production 
(decay) operator 

J{x) = y e-^*^^ b.V^ci W^ciXc2] {x) + h.c. (50) 

Here we have added the Wilson fine (see (47)) to the deflnition of the operator. This 
generates an infinite number of interaction vertices with additional n+Ad photons fields, 
which are not suppressed by powers of 6, but all have the same short- distance coefficient C 
because of gauge invariance [15]. The Wilson line accounts for full-theory diagrams where 
a soft 0-field emits one or several cl-coUinear photons and, thereby, becomes off-shell. 
Integrating out the hard internal 0-field results in the Wilson fine. At next-to-leading 
order we are only concerned with the emission of at most one coUinear photon, which must 
then couple to the external electron line attaching to the same production (decay) vertex 
from which the coUinear photon is emitted. We also mention that in principle we should 
consider the emission of photons coUinear with the neutrino from the electron or the heavy 
scalar, which puts the electron or scalar Unes off-shell. However, the effect cancels due 
to the charge- neutrality of the neutrino. It is for this reason that we do not introduce a 
c2-collinear photon field in the effective theory. 
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We begin the one-loop computation of C by specifying the matching prescription. In 

the underlying theory we compute the amputated, MS-renormalized, on-shell three-point 
function of an electron, neutrino and scalar field, multiplied by the LSZ residue factors, and 
expand it to leading order in 5. The corresponding three-point function in the effective 
theory contains one insertion of the operator y e~'^^'"^ [(t)vipciWciXc2]{x). Denoting this 
three-point function by Feff and the former by F, the matching equation reads 



Here the various i?-factors denote the ultraviolet-finite LSZ residues in the MS scheme 
given by the ratios of the on-shell and MS field renormalization constants. The factor 
defined in (32) appears because of the non-relativistic normalization of the field 0^. 
Its origin can be understood most directly by comparing the scalar propagators in full 
and effective theory diagrams as done in (31). Note that the on-shell condition for the 
scalar line implies that its momentum satisfies = s = -|- MA, but in a perturbative 
matching calculation this condition must be fulfilled only to the appropriate order in a 
and 5. For the one-loop matching of C it is sufficient to put = M^. 

The matching equation (51) is valid for any infrared regulator, but becomes particularly 
simple in the case of dimensional regularization, where all loop diagrams in the effective 
theory vanish on-shell. Then all the effective theory residues equal 1, Fgff equals its tree- 
level value and the matching equation becomes 



The full theory quantities must then be equal to their hard contribution in the sense of 
the expansion by regions, and since w can be expressed in terms of the short-distance 
coefficient A, this equation says that we can obtain C by directly computing the hard 
contributions in the full theory. 

The residue Rhrj, can be obtained from (35) and (39), and the residues P/i^ and Rhx 
can be computed in a similar way. The results are given in the Appendix. There is only 
one one-loop diagram and one counterterm diagram contributing to F at next-to-leading 
order. Evaluating the hard part of this one-loop diagram we obtain 




(51) 




(52) 




(53) 




(54) 



and the residues from the Appendix we obtain C — C^^^ -\- O {o?), where 
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Figure 3: Hard contributions to the next-to- leading order amplitude. Left: Inser- 
tion of [A(^^]^/4 and — MA^^^. Right: Insertion of C^^\ the symmetric diagram is 
understood. 



The poles are infrared poles related to the factorization of hard and soft contributions. The 
final result for C^^^ is obtained by subtracting the pole part. (In the following C'-^^ refers 
to this subtracted expression.) It is independent of the gauge parameter ^ as it should be, 
because the production (decay) operator is gauge- invariant. The double logarithm in the 
coefficient of ag is related to the infrared divergence in F^. The single logarithms on the 
other hand are also related to the scale dependence of the Yukawa coupling y. 



3.4 Forward scattering amplitude in the effective theory 

We are now in the position to compute the forward scattering amplitude at next-to-leading 
order in the effective theory. 

The first class of contributions consists of diagrams with the topology of the tree di- 
agram (22), but with corrections to the propagator and the production (decay) vertex as 
shown in Figure 3. As can be seen from (44), the relevant terms for the left diagram 
of Figure 3 are [A(^^]^/4 and —MA^'^\ The remaining term (zDst)^ does not lead to a 
contribution at NLO. First for head-on collisions the scalar line has kj = 0. Secondly, the 
photon vertices from this term contribute only at NNLO for weak coupling g. The other 
diagrams in this class are obtained from the hard corrections to the Wilson coefficient 
C of the production (decay) vertex. (The production and decay operator have the same 
complex coefficient C, and not the complex conjugates to each other, see (50).) There are 
two such diagrams, one of which is displayed in Figure 3. Adding up these three diagrams 
we obtain 



8VM 2V 



(56) 



with 

. aw 

V = ^-M —. (57) 

All these diagrams correspond to corrections due to hard loops and, thus, can be identified 
with the factorizable corrections [10]. These corrections are separately gauge invariant 
since so are the matching coefficients. In addition there is the contribution from the 
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Figure 4: Soft contributions to the next-to-leading order amplitude. A diagram 
with a soft correction at the decay vertex is understood. 

production-decay operator T in (49), which contributes to iT^^^ the term 




iyy_s_ 
8M2 
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We now turn to the one-loop corrections in the effective theory, which are due to soft 
or coHinear photons. The loop provides a factor of a, so at NLO in a combined expansion 
in a and 6 we can use the leading power soft and coUinear photon vertices. We begin with 
the soft correction. The coupling of the soft photon to the 0^-field is given by the covariant 
derivative in the first term in £hset (44), whereas the coupling to the coUinear electron 
can be obtained from the first term in £scet (7). With dimensional infrared regularization 
the residue factors equal 1, which leaves the diagrams shown in Figure 4. For the forward 
scattering amplitude the box diagram (third diagram in Figure 4) vanishes in Feynman 
gauge, because the photon-coupling to the coUinear electron is proportional to n'^. Thus 
the ampUtude is proportional to =0. Computing the remaining diagrams of Figure 4 
we obtain 



X a„ 



-2e 



4 + 



6 



(59) 



The final result for the soft amplitude is the above expression with the pole terms subtracted 
after expanding the prefactor (— 2D//i)"^^ The soft part of the amplitude does not contain 
logarithms of the hard scale as expected. We also note that the soft amplitude is gauge- 
invariant, because it is computed with a gauge-invariant effective Lagrangian. 

The only coUinear photon correction at NLO comes from the diagram shown in Figure 5 
(and the corresponding symmetric diagram), where the coUinear photon coupling at the 
production (decay) vertex is derived from the Wilson line in (50). However, for an on-shell 
electron the integral is scaleless and vanishes in dimensional regularization. 

Combining (56), (58), (59) with the result for the matching coefficient C^^-* in (55) we 
obtain for the forward scattering amplitude at next-to-leading order 
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Figure 5: Contributions to tiie next-to-leading order amplitude involving collinear 
photons. 

/ -2MV 3, M2 7 27r2\ 

+ 2 log — - - - ?7r * 



M /i^ 2 / 8PM 2D 2M 



(60) 



Rather than inserting the subtracted expressions for T^*-^-* and C^^-* we have used here the 
unsubtracted ones to display the structure of singularities. The double poles in 1/e have 
canceled between the hard and soft corrections, but there are 1/e poles left over from T^^\ 
which must be associated with the collinear singularity due to initial state radiation from 
the massless electron. To verify that these singularities have the expected structure we 
compute 

[^^sS2/^] = Cdx-P{x) [zrW(a:s)/(x.)] (61) 

ZTT JO e 

where iT^^\xs) is the leading-order forward scattering amplitude given in (22), but with s 
replaced by xs. {x is the momentum fraction of the electron in the electron.) Furthermore 

is the soft limit of the e — > e Altarelli-Parisi splitting function, and we took into account 
that the incoming neutrino does not radiate photons. Simplifying the splitting function 
to the soft limit is justified, because values of x not near one correspond to hard-coUinear 
radiation, and this is a NNLO effect (which is part of the production-decay operator in the 
right panel of Figure 1). Similarly the x-integral could be restricted to an interval close 
to one. Equivalently, we perform the x-integration from zero to one and expand the result 
around ^/s = M to obtain 

3 + 41ogZ?a£5V (63) 

M2 y 

which agrees with the pole part of (60). This confirms that the cross section is finite once 
the initial-state singularities have been factored into the electron distribution function. 
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Subtracting the initial-state singularities minimally according to our convention for the 
electron distribution function, we obtain 



X 



/ , -2MV , -2MV ^ -2MV 
ag 3 log h 4 log . log 



-2MV 3, 7 27r2^ 

' ^og~rv; o log — ~ o + 



M2 2 ° 2 3 y 

+ I 21og — - - - ZTT - '-—4- + 




(64) 

2 ; 8PM 2V 2M ^ ' 

as the final result for the next-to-leading order contribution to the forward scattering 
amplitude. The correction to the inclusive line shape is given by the imaginary part of 
this expression. We have introduced the scale v to distinguish the scale of the electron 
distribution function from the renormalization scale in the masses and couplings. Referring 
to (25) we note that v should be taken of order M5 to make the initial-state coUinear 
logarithms small. 

We performed several further checks of this result. First we find that the logarithms 
of 11 cancel the scale dependence of M and yy* in iT^^^ to NLO. Second we expand 
1/s X Im (T*^°) -|- T*^^^) to order cx^ and compare this to the expansion of the fixed-order 
cross section (25) to order The two expansions coincide in the terms o? jP'^ 

which they have in common. Finally, we have also computed the 5 expansion of the one- 
loop forward scattering amplitude directly in the full theory. This checks all the terms in 
(64) except for those involving A*^^) ^^(1 [A*^^)]^, which come from two- loop diagrams in the 
full theory. 

The final result (64) is given in the MS scheme. We briefly discuss how to translate this 
result to the pole renormalization scheme, where the mass M is defined by s = - iMV 
and the residue of the renormalizcd (^-propagator is i?^ = 1. All other renormalization 
conventions remain unchanged. From the definition of A in (28) we obtain in the pole 
scheme Re Ap = and 

Im^^ Im(A/M) 

M l + Re(A/M) 

with A and M referring to the MS (or, in fact, any other) scheme. Using (42), (43) we 
obtain 

A(i) 

-jf = -2n:ay, (66) 

= ^7ra2(-61n— + 5j+z7ra,aJ61n— -27+^j. (67) 

To convert the matching coefficient C to the pole scheme, we observe that R^Z^j, is scheme- 
independent, so the only modification comes from the normalization factor ro, resulting 
in 

C,= {^\'"C, 2CW = 2C« + 5^. (68) 
^ \wj ' ^ 2M 
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This is exactly what is needed to render the hne shape scheme-independent to leading 
order in S and to all orders in a. Using (65) and (68) together with — + MRe A 
we find C^/M ^C^M, so 

= (69) 

exactly. 



3.5 Next-to-leading order line shape 

We compute the line shape at NLO in the MS renormalization scheme by taking the 
imaginary part of (T*^°^ + T^^^)/s. The result depends on the renormalization scale and 
rather strongly on the collinear factorization scale u. This dependence would be canceled 
(up to the order at which we perform the calculation), if we folded our cross section with 
the electron distribution function. Since this is well understood and we are interested 
in displaying the next-to-leading order effects intrinsic to the line-shape calculation, we 
present the partonic cross section, and adopt fi = M and u = \2V\ unless otherwise stated, 
since this choice makes the collinear logarithms unimportant. 

In order to illustrate the numerics we use Qg = 0.1/(47r), ay = 0.1 /(47r) and ax — 
(0.1)^/(47r)^ and assume the pole mass M — 100 GeV. For the plots in the MS scheme we 
have to convert the pole mass to the MS mass M. At LO we use the one-loop relation 
between the pole mass and M and obtain M = 98.8 GeV, whereas at NLO we use the 
two-loop relation and get M = 99.1 GeV. The LO line shape has already been shown in 
Figure 2. In Figure 6 (upper panel) we show the LO and the NLO hne shape in the MS 
scheme (solid curves) and the pole scheme (dashed curves). The scheme dependence is 
very small, but the next-to-leading order correction is around —10% near the peak, and 
can be up to —30% for center- of- mass energies between 95 GeV and 105 GeV . The ratio 
of NLO to LO is shown as the solid line in the lower panel of the figure. 

Contrary to the LO hne-shape the NLO line shape does not any longer have an exact 
(non-relativistic) Breit-Wigner shape. Performing a fit of the NLO hne shape (MS scheme) 
to a Breit-Wigner curve in the range 95 < ^/s < 105 GeV we find deviations up to 
15%. The ratio of the NLO curve to the Breit-Wigner fit is shown in as dashed line 
in Figure 6. Moreover, the output mass parameter of the Breit-Wigner fit differs from 
the input M = 100 GeV by — 160MeV. This illustrates that in realistic situations the 
data should be fitted to the theoretically predicted shape rather than a Breit-Wigner-type 
ansatz. 

Up to now we have performed a strict expansion in the couplings and S to obtain the 
NLO result. In particular, we only included A*^^^ in the propagator and treated the terms 
A*^^) as well as higher powers of A*^^) as interactions in a perturbative expansion. In other 
words, we used i/{2MT>) for the propagator of the unstable field. A shghtly different 
approach is to include the effect of higher-order bilinear terms in the Lagrangian into the 
propagator, which sums these terms to all orders. This can be done in different versions. 
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Figure 6: Upper panel: line shape (in GeV~^) in the MS scheme (solid) and 
pole scheme (dashed) at LO (upper magenta/light grey curves) and NLO (lower 
blue/dark grey curves) as a function of the center-of-mass energy in GeV. Lower 
panel: The ratio of the NLO to the LO line shape in the MS scheme (solid blue/dark 
grey curve) and the ratio of the NLO MS line shape to a Breit-Wigner fit (dashed 
magenta/light grey curve). 

For example, instead of A*^^) we may use the full A. At NLO, this implies using the 
propagator 



2M(P-A(2)/2) ^'"^ 

which absorbs the interaction term bS'^^ jTD and sums all powers of A*^^-* in the forward 
scattering amplitude. The difference to the non-summed line shape is of course a NNLO ef- 
fect. The ratio of summed to non-summed line shape at NLO is shown as the dashed curves 
in Figure 7 in the MS scheme (blue/dark grey line) and the pole scheme (magenta/light 
grey line). In computing the resummed line shape we have adapted the factorization scale 
and used v = \2V - A(2)|. It can be seen from Figure 7 that the effect of resummation 
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Figure 7: The ratio of the line shape with powers of A^^^ (and A*-^^) resummed over 
the non-summed line shape in the pole (magenta/light grey lines) and MS scheme 
(blue/dark grey lines). The dashed curves represent the line shape including A^^^ 
in the propagator. Resumming also all powers of A^^^ results in the solid curves. 

is less than 2%, which indicates that the higher-order terms associated with A*^^-* are not 
large. The effect is still smaller in the pole scheme. 

We can go one step further and also resum all terms associated with the expansion of A. 
In order to achieve this we have to recall that we obtained the propagator by expanding 

vk = —M + \J + MA — kj in A and kj (see (29)). Expanding only k^ we obtain the 
summed propagator, which at NLO reads 



2M( - ^JM^ +M{AW + A(^y)) ' 

This also absorbs the interaction term — [A*^^)]^/8MD into the propagator. The ratio of 
corresponding summed to non-summed line shape at NLO is shown as the solid curves in 
Figure 7 in the MS scheme (blue/dark grey line) and the pole scheme (magenta/light grey 
line). (The factorization scale has been adapted to z/ = 2\{^/s-yjM'^ + M(A(i) + A(2)))|.) 
Again the effect of resummation is marginal compared to the size of the NLO correction. 

4 Beyond NLO 

A crucial aspect of our formalism is that it allows to improve the accuracy of the calculations 
by including systematically higher-order corrections. After the explicit calculation of the 
line shape at NLO, we now discuss how to go beyond NLO. We will outline the necessary 
calculations to obtain the line shape at NNLO in the effective theory approach. Then 
we compute explicitly a subset of corrections to illustrate once more how the formalism 
automatically maintains gauge invariance of the calculation. 
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4.1 Elements of the NNLO calculation 

Since we perform a double expansion in a and 5, the NNLO computation must include all 
terms suppressed with respect to the LO result by either a^, 5^ or a5. In order to facilitate 
the classification of terms we also distingTiish between radiative corrections due to hard (h), 
collinear (c), and soft (s) loops, denoting the coupling a with the corresponding subscript. 
For example, a two-loop matrix element correction in the effective theory with no further 
suppression by 5 will be denoted as a^, if both loops are soft, and if the vertices do not 
include hard corrections. With this in mind we proceed to the discussion of various NNLO 
contributions. 

1) Tree amplitude in the effective theory with insertions of LO operators in 5, matched to 
NNLO (al) 

The effective theory diagrams are the same as in Figure 3, but now the propagator and hard 
vertex corrections must be included to higher accuracy in a. This implies the computation 
of A(=^), which, from (28) and (34), can be seen to involve n(^'°\ the three-loop self-energy 
at = M^, as the most difficult part. In addition, the two- loop hard correction C^^^ to 
the production (decay) vertex (j)v'4'ciXc2 is required. These are standard loop calculations, 
which arc difficult, but not impossible with present Feynman diagram technologies. Both 
matching coefficients, A^^^ and C*^^-*, are gauge invariant. 

2) Tree amplitude in the effective theory with insertions of NLO operators in S, matched 
to NLO (ahS) 

The number of possible terms of this type is limited by the requirement that they con- 
tribute to the tree 2 — > 2 scattering amplitude in the effective theory. NLO operators 
can be constructed from the LO operators by acting on soft and collinear fields with extra 
derivatives. The Wilson coefficients of these operators then have to be matched up to NLO. 
Operators with bilinear terms in the fields require no work, since they are determined by 
the dispersion relation for the free field. Operators generated by adding derivatives to the 
LO production (decay) operator can be eliminated by the equation of motions or do not 
contribute at tree level for external states with vanishing transverse momentum. 

This leaves genuine NLO operators. An example is the four-fermion production-decay 
operator (49), 

(4lXc2)(Xc2V'cl), 

4M2 

that we encountered already in the computation of the NLO scattering amplitude. At 

NNLO we therefore need the hard radiative correction to the coefficient function. The 
corresponding calculation is done explicitly in the subsequent subsection. In general, pro- 
duction (decay) operators need not be local. Consider a collinear electron fluctuation with 
momentum Mn_/2 -|- ki (where n+ki is of order M for a collinear fluctuation) scattering 
on the neutrino with momentum Mn+/2 -|- k2, where k2 is small. In this case, the scalar 
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propagator is i/{Mn+ki), which is far off-shell but does not correspond to a local inter- 
action vertex. However, as discussed before, these non-local operators do not have tree 
matrix elements with the external collinear modes in 2 — > 2 scattering. 

Finally, there is a one-loop correction to the four-fermion production-decay operator 
that arises indirectly. Loop corrections induce an effective soft photon-neutrino coupling 

^ (Xc27.Xc2) d,Fr (72) 

Although this operator scales as gayS^, which is gayS^ suppressed relative to the neutrino 
kinetic term, it contributes to the NNLO forward scattering amplitude through soft-photon 
exchange in the ^-channel, because the soft photon propagator cancels the 5^ suppression. 
This can be seen directly by converting (72) to the four-fermion operator using the soft 
photon equation of motion. In the following subsection we compute the coefficient function 
and discuss how this operator contributes to the tree scattering amplitude in the effective 
theory. There exists also a similar electron-photon coupling. However, in this case there is 
no tree scattering diagram due to the absence of a tree-level neutrino-photon coupling. 



3) Tree amplitude in the effective theory with insertions of NNLO operators in 5, matched 
to LO (5^) 

The task of writing down all relevant 5"^ suppressed operators is simplified by the fact that 
we do not consider further suppressions by a. Hence we only need to expand the tree 
2—^2 scattering diagram in the underlying theory one order further than in (48). We 
then find the NNLO production (decay) vertex '?/^ciXc2(^-Dst)^0?; (which in fact gives no 
contribution in the center-of-mass frame of the coUision, where k-j — Qi), and the NNLO 
production-decay operator 

i'^clXc2) ivDs (Xc2'0cl), (73) 

which produces a correction proportional to (y^ — M)/M^ to the line shape. 



4) One-loop amplitude in the effective theory with insertions of LO operators in 5, matched 
to NLO (agOth^oicOih) 

The one-loop scattering diagrams are those shown in Figure 4, however, now there is an 
additional insertion of [A(^)]^/4 or — MA^^^ into the scalar propagator, or the production 
(decay) vertex is taken at NLO, involving the correction C'^^\ similar to the tree diagrams in 
Figure 3. These terms can be obtained from the calculation in the previous section without 
further work. We simply muhiply (59) by interpret V as ^-(M2 + M(A(^) + A(2)))^/2^ 
expand in a and 5 and pick out the NNLO terms. Note that in the underlying theory 
the two-loop diagrams corresponding to these contributions involve soft and hard loops, 
hence in the conventional terminology these are neither purely factorizable nor purely 
non-factorizable corrections. Nevertheless, the effective theory formalism guarantees that 
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Figure 8: Contribution to the next-to-next-to-leading order amplitude involving 
collinear (and soft) photons. 

this set of terms is gauge-independent. Since the colhnear one- loop contribution shown in 
Figure 5 vanishes, it follows that there are no terms proportional to acUh- 

5) One-loop amplitude in the effective theory with insertions of NLO operators in 5, matched 
to LO (ois5, ac5) 

We have already discussed possible NLO operators in the category 2), but some of the 
simplifications based on the properties of the external momenta of the tree amplitude in 
the effective theory no longer apply to the one-loop amplitude. For instance, the residual 
momentum fcj of the internal scalar line no longer vanishes, and a non-zero contribution 
from the insertion of {iDg-jY into the scalar propagator is obtained. In principle one may 
also need new production (decay) operators with collinear or soft photon fields in addition 
to the basic structure (f)v'4'ciXc2- The problem of deriving the operators at tree- level is 
very similar to deriving the power-suppressed heavy-to-light decay current in soft-collinear 
effective theory [16]. Since the collinear one-loop diagrams vanish, the interesting operators 
are those with a soft photon field, but it turns out that after applying the field equations, 
there is no gauge- invariant operator of this sort. It follows that the set of possible NNLO 
terms is again rather limited. One as6 term that appears is the soft one-loop correction to 
the local four-fermion interaction shown in (58). 

6) Two-loop amplitude in the effective theory with insertions of LO operators in 6, matched 
to LO (a1, asexual) 

The final set of terms involves a two-loop calculation of the forward scattering amplitude 
in the effective theory, including vertex- and box-type corrections. However, since the 
propagators and vertices of the leading-order effective Lagrangian are much simpler than 
those in the underlying theory, the two-loop calculation is also greatly simplified. It is 
worth noting that in this calculation there appear for the first time non-vanishing collinear 
loop integrals. A representative diagram is shown in Figure 8. The requirement is that the 
external electron first emits a soft photon, so that the subsequent collinear emission occurs 
from an off-shell line. The collinear loop integral is then no longer scaleless. Because of 
this requirement, however, there are no contributions where both loops are collinear. 
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This concludes the classification of terms. It is evident that a NNLO calculation of 
the line shape can be performed in the effective field theory approach without conceptual 
difficulties. In fact, the previous discussion shows that the most difficult calculations are 
reduced to standard calculations of a hard two-loop vertex correction and the on-shell three- 
loop self-energy. In particular, the difficulty of maintaining gauge invariance is absent in 
this approach. Since the separation of terms by momentum scales is gauge- independent, 
the sum of all contributions with a given power of 5, ah, ckc, and cts is separately gauge- 
independent. 

In the version of the effective theory where coUinear fluctuations are integrated out, the 
coUinear loop contributions move from the effective theory matrix element to a matching 
correction. The calculations to be performed remain of course the same, but the structure 
of operators is different, since the effective theory no longer contains fields for coUinear 
fluctuations. While this simplifles the fleld content of possible operators, it also complicates 
their structure, introducing new non-localities related to the fact that one component of 
coUinear momentum is of the same order as the soft momentum. There is again a close 
analogy with recent developments in soft-coUinear effective theory, which, however, we do 
not pursue further here. 

4.2 One- loop matching of the production-decay operator 

In this subsection we perform an explicit calculation of the ahS contributions to the forward 
scattering amplitude categorized under 2) above. 

4.2.1 Four-fermion operator 

The four-fermion operator (49) is a sub-leading operator that gives rise to a production- 
decay vertex. It contributes to the forward scattering amplitude at NLO through the tree 
diagram (58). We now determine the loop correction to the coefficient function of this 
vertex. 

Consider the box diagram contribution to the scattering amplitude in the full theory 
shown in Figure 9. It is clear that the hard part of this diagram results in a loop correction 
to a local four-fermion operator. This contribution alone is gauge-dependent. However, 
there are additional contributions from vertex and self-energy corrections (also shown in 
Figure 9), since expanding the corresponding subgraphs in 5 results in terms where the 
propagators of the unstable scalar are canceled. Because of this cancellation these terms 
also contribute to the coefficient of the local four-fermion operator, rendering the complete 
result gauge-independent. 

To see this, consider first the terms proportional to ayOg in Figure 9. Due to the 
two scalar propagators, the leading contribution from the self-energy diagram is of order 
o? jP. This is a LO contribution, which is already included in the effective theory through 
the propagator i/{2M{yk) — MA*^^-') of the scalar field. The next term in the expansion 
in 6 cancels one propagator, resulting in a local vertex. This is already included in the 
computation of the matching coefficient C'-^-*. The second term in the expansion in 5, 
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Figure 9: Diagrams in the full theory proportional to OyOg, whose hard part con- 
tributes to the matching coefficient of the production-decay vertex at NLO. A dia- 
gram with a vertex correction at the decay vertex is understood. 

however, results in a NNLO contribution, which has to be reproduced by the effective 
theory via the matching coefficient of the production-decay vertex. Similarly, the leading 
contribution from the vertex diagram is included in C^^\ and we need the next term 
in the expansion in 6, where the scalar propagator is canceled. We then obtain for the 
contribution from the self-energy, vertex and box diagrams 



^(2) _ . ,2 Hp)v{q)] [v{q)u{p)] .^(1,2) 



'-^Hp)vmmu{p)] + + (74) 



z|yPa„r-/ N / M r-/ N / M / 2 2-^ 2, M2 

Hp)v{q)] [v{q)u{p)\ — - + + - In — 

\ ^ e e jj. 



M2 , 2M2 71^' 



(2-01n^-ln^^-2(4-0-y , (75) 



Sf = ^[«(p)^(9)][%Mp)] (-^ + ^ln^ + (l-0 



Summing (74), (75) and (76) we see that the gauge dependence cancels, and we obtain 

^-S^[u{p)v{q)] [v{q)u{p)] + ^-^[u{p)-f''Yviq)] [v{qhu-ff,u{p)] (77) 



with 



, ,2 /2 2, M2 , 2^^ 7r2\ 
5r = -M% ---ln- + ln^- + 6 + - , 



, ,2 / 1 1, M2 3\ 
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Figure 10: Diagrams in the full theory proportional to a^, whose hard part con- 
tributes to the matching coefficient of the production-decay vertex at NLO. 



The cancellation of the gauge dependence is of course not accidental. Since the separation 
of hard and soft parts is gauge- invariant, and since we match S'-matrix elements, the 
short- distance coefficient must be gauge- independent. 

There are further contributions to the matching coefficient. First, the terms propor- 
tional to ay from the diagrams shown in Figure 10 are given by 

iB iB iB 

-^[uiPHq)] [v{q)u{p)] + -^[u{p)(iu{p)] [v{q)M(l)] + -^[u{p)Yu{p)] [v{q)-f^v{q)] (79) 
with 

Bs = -\yW, B, = -\yW (^-5 + , B, = -\yW (l " ^) ■ (§0) 

Second, up to now we have been discussing contributions which (with the exception of 
the box diagrams) come from expanding the one-particle irreducible self-energy or vertex 
subgraphs in 6. Additional contributions arise from the expansion of the one-particle 
reducible scalar propagator according to (48). For instance, combining the local NLO 
term from (48) with the (renormalized) hard vertex at leading order in 6, we obtain a 
contribution to the coefficient of the four-fermion operator proportional to C^^\ Similar 
terms come from the self-energy diagram. Collecting all contributions, we obtain for the 
NLO correction to the production-decay operator (49) 



T^'^ = -^(^clXc2)(Xc2^^cl) 



a, + aJ---ln— + ln ^ + y + 6 



—r- + 



16M 



y a u , ^ / 1 1 3 

VPcilvlt,Xc2) {Xc2Tl ^ci) a J — ~ o In — + 



Xc2 ) ay 




\y\ - - / ^ 

{ipcll^lpcl) {Xc2l^Xc2) S I 1 - 
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Figure 11: NNLO forward scattering diagrams in the full theory with a soft pho- 
ton exchanged in the t-channel. The hard contribution to the triangle subgraph 
generates the effective photon-neutrino coupling {xx)dF. 



with A(i)/M and given in (42) and (55), respectively. Inserting the explicit expres- 
sions for these quantities, the poles in e do not cancel. The left-over poles are related to 
the factorization of hard, coUinear and soft contributions, and the initial-state collinear 
singularity, and are canceled only once all other NNLO contributions to the scattering 
amplitude are combined. The final result for the matching coefficient is obtained from the 
above expression with all poles subtracted minimally. 

We expressed (81) in terms of four-fermion operators with four Dirac structures as they 
come out of the calculation. Not all of them are independent. Using Fierz transformations 
and the projection properties J^-V'ci = ^+Xc2 = of the collinear fermion fields, they can 
be reduced to linear combinations of the three basis operators 

ii>clXc2)iXc2^cl), (^cl75Xc2)(Xc275^cl), (^c17a<x Xc2) (Xc27^^ ^cl) , (82) 

where "/ij," means that one sums only over the transverse components. To perform this 
reduction, it is necessary to define carefully the Dirac algebra in d dimensions in very much 
the same way as in applications of the weak effective Hamiltonian to weak interaction 
processes at low energies. 



4.2.2 Effective photon-neutrino coupling 

There is a second correction of order atS, which we discuss separately, although it can be 
interpreted as a further contribution to (81). Consider the two diagrams in the full theory 
displayed in Figure 11, which do not involve resonant scalar lines. The hard contribution 
to the triangle subgraphs (consisting of scalar and electron lines) induces an effective soft- 
photon neutrino coupling. The result of the calculation gives 



S (Xc27.Xc2) d,Fr 



1 

3^ 



1 M2 1 

- In — + - 

3 ^2 2 



^3) 



In Figure 11 the photon must be soft, because we consider the forward scattering amplitude. 
(In fact, to regulate the otherwise divergent t-channel propagator, we must temporarily 
consider the amplitude slightly non- forward.) Power counting shows that this is a 5^ 
suppressed operator. The factor of 6"^ is compensated by the t-channel photon propagator 
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Figure 12: Effective tlieory tree diagrams related to tlie full theory diagrams in 
Figure 11. The left-hand diagram involves the {xx)dF coupling, the right-hand 
diagram the mixed collinear-soft four-fermion operator and a soft electron loop. 

of order 1/6^, the result being an agOy (NNLO) contribution to the scattering amplitude in 
the effective theory shown as the left-hand diagram of Figure 12. The photon propagator 
is canceled by a momentum factor from the vertex (83) resulting in a local contribution to 
the scattering amplitude of the same form as (81). This can be seen explicitly using the 
soft photon equation of motion d^F^!-" = —g ipcil^i^ci + ... to rewrite (83) as a four-fermion 
operator times a coefficient proportional to \y\'^ag. 

The 1/e pole in (83) is a soft divergence, which arises due to the hard-soft factorization 
from the second diagram in Figure 11. The divergence cancels with the soft contribution 
to this diagram, where the two internal electron lines are soft, and the scalar line is off- 
shell. The off-shell scalar induces a local {Xc2'ips){'4'sXc2) vertex. This mixed collinear- 
soft four-fermion operator contributes to the forward scattering amplitude in the effective 
theory through the one-loop diagram shown on the right-hand side of Figure 12, where the 
electron loop is soft. Once again the 6^ suppression of the operator is compensated by the 
soft photon propagator. 

5 Summary and outlook 

It has been recognized for some time that the perturbative treatment of unstable particles 
is difficult, because partial summations of the perturbative expansion are necessary. These 
difficulties originated mainly from the fact that the guiding principle for resummation was 
not understood. In this paper we advocated the idea that effective field theory methods 
familiar from other applications in high-energy physics can be adapted to solve the prob- 
lem of describing systematically the production and decay of unstable particles, since the 
fundamental reason for the breakdown of weak-coupling perturbation theory is related to 
the emergence of a second (small) momentum scale near resonance. The advantages of the 
effective theory method are: 

• It breaks the calculation into several well-defined pieces (matching calculations, ma- 
trix element calculations), thus rendering the organization of the calculation efficient 
and transparent. 

• It provides a power-counting scheme in the small parameters {S, couplings), which 
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allows for an identification of the terms relevant for achieving a prescribed accuracy 
before actual calculations must be done. 

• It provides a set of (Feynman) rules to compute the minimal set of terms necessary 
for a given accuracy. Since one does not calculate "too much" , the calculation to a 
given order is presumably technically simpler than in any other approach. 

• Gauge invariance is automatic at every order, since the effective Lagrangian is gauge- 
invariant. 

• It can be extended to any accuracy in the expansion in S and in couplings at the 
expense of performing more complicated, but well-defined calculations. 

In this paper we explained these ideas in a toy theory restricting ourselves to a single scalar 
resonance, an abelian gauge theory, and a totally inclusive resonant scattering cross section. 
Extending the method to non-abelian gauge theories or relaxing any of the other restric- 
tions causes only technical complications. For instance, for a resonant fcrmion, the heavy 
scalar Lagrangian would be replaced by the familiar heavy quark effective Lagrangian. For 
resonant gauge bosons in the clcctrowcak Standard Model in gauge with ^ not near 
one, the degrees of freedom with mass are integrated out, and the effective theory 
contains a massive vector field with three polarizations. Non-renormalizability is not an 
issue, since the effective field describes only soft fiuctuations. In 't Hooft-Feynman gauge, 
or with ^ close to one, the efi^ectivc theory contains a vector field with four polarizations 
and a pseudo-Goldstone field, which cancels the effect of the scalar polarization state, so 
that one can work again with a canonical massive vector field. 

Other obvious extensions of this work concern non-inclusive kinematics and pair pro- 
duction of unstable particles. In these cases it will be necessary to enlarge the field content 
of the effective theory. For instance, if energetic particles (jets) are detected in the fi- 
nal state, one must introduce new collinear fields corresponding to the direction of these 
particles. It is advantageous to work with cut diagrams rather than amplitudes, so that 
the (final-state weighted) phase-space integrals can be treated on the same footing as loop 
diagrams in matching calculations. Similarly the extension to pair production requires the 
introduction of two copies of what we called the HSET Lagrangian in this paper, one for 
each unstable particle with different velocity vectors. The production (decay) vertices then 
depend on additional kinematic quantities such as the scalar product of the velocities. Pair 
production near threshold appears in this context as the particular case, where the two ve- 
locities become nearly equal. The two HSET Lagrangian then merge to the non-relativistic 
Lagrangian, and standard methods can be applied to deal with the potential complications 
due to the strong Coulomb force. We are confident that effective field theory will be 
the method of choice to address unstable particle production in these different kinematic 
situations and plan to return to some of them in the context of concrete applications. 
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Appendix 

In this appendix we collect the renormalization factors and counterterms of the theory 
defined by the Lagrangian (2). The counterterm Lagrangian is 

+ 5yyi,^cl>^l;x + S;y*i^'<l>^xi^ - (84) 
The relation between bare and renormalized quantities is 



00 = 






Z<t, — 1, 






5m = 


- Zl^Zfj, — 1, 


V'o = 






Z^j, — 1, 


Xo = 




Sx = 


^x-l = 


Ao = 




5a ^ 


Za-1. 


^0 = 








9o = 


Zggfi'', 


^ = 


ZgZ^\J Za — 1, 


yo = 




5y = 


Zy\J~Z^^J~Zp^J~Z^ — 1, 


Ao = 


(A + A,)/.^^ 


Sx = 


(Z|-1)A + Z|A,, 



(85) 



The counterterms 5'g and 5g2 are not independent, but determined by the previous ones: 



5g — ZgZ^\j Za — 1, 

b'g, = ZIZ^Za-1. (86) 

In this paper all renormalization constants are needed at the one-loop order except for 
which must be known to two loops. In the MS scheme, and in covariant gauge with gauge 
parameter the counterterms expressed in terms of the renormalized couplings are given 
by: 

1/3 3 

Zy = 1 + - (^-2«9+ g"^' 
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= l + ^((3-Oa,-2a,) 

Z^ = 1 + - (-^(Ig -^dy 



Za = ! + -(--% 



e V ^ 2 
e V 2 

-f-- 

e V 3 



Zl^ = 1 + - (-3 + 2 Oy) + ^ (S + - 9 a^^a^^) 
1 /^53 2 3 2 5 \ 



Furthermore, we need the part of ?s\ that is proportional to o?y and a^. It is given by 



e e 

The ultraviolet-finite residues of the MS-renormalized propagators read: 

1 / \ 

= l-(3-0-a3 + «y |21og— -2-2i7rl 



= 1 + 02/ (-log 



1 M2 r 



/ 1 4 \ 

= l + „,(^_l„g_ + _j (89) 

The 1/e poles in these equations are infrared divergences. In the dimensional infrared 
regularization scheme the full residues equal their hard contributions, so R^x — Rx- 
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